Conditions for free magnetic monopoles in nanoscale square arrays of dipolar spin ice 
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We study a modified frustrated dipolar array recently proposed by Moller and Moessner [Phys. 
Rev. Lett. 96, 237202 (2006)], which is based on an array manufactured lithographically by Wang et 
al. [Nature (London) 439, 303 (2006)] and consists of introducing a height offset h between islands 
(dipoles) pointing along the two different lattice directions. The ground-states and excitations are 
studied as a function of h. We have found, in qualitative agreement with the results of Moller and 
Moessner, that the ground-state changes for h > hi, where hi = 0.444a (a is the lattice parameter or 
distance between islands). In addition, the excitations above the ground-state behave like magnetic 
poles but confined by a string, whose tension decreases as h increases, in such a way that for h « hi 
its value is around 20 times smaller than that for h — 0. The system exhibits an anisotropy in the 
sense that the string tension and magnetic charge depends significantly on the directions in which 
the monopoles are separated. In turn, the intensity of the magnetic charge abruptly changes when 
the monopoles are separated along the direction of the longest axis of the islands. Such a gap is 
attributed to the transition from the anti to the ferromagnetic ground-state when h — hi. 

PACS numbers: 75.75.-c, 75.60.Ch, 75.60.Jk 



INTRODUCTION 



Geometrical frustration in magnetic materials occurs 
when the spins are constrained by geometry in such a 
way that the pairwise interaction energy cannot be si- 
multaneously minimized for all constituents. A special 
example is an exotic class of crystalline solid known as 
spin ice (Dy2Ti20?, H02T12O7). Recently, Castelnovo 
et alH have proposed that these materials are the repos- 
itory of some elegant physical phenomena: for instance, 
collective excitations above its frustrated ground-state 
surprisingly behave as point-like objects that are the con- 
densed matter analogues of magnetic monopoles. Some 
recent experimentsEra have reported the observation and 
even the measurement of the magnetic charge and cur- 
rent of these monopoles in spin ice materials; in addition, 
simulations also support these ideasOil Besides, to turn 
the research of monopoles into a proper applied science, 
it will be necessary to ask if the basic ideas of dipole 
fractionalizationBH that give an usual spin-ice material 
its special properties can be realized in other magnetic 
settings. One of the most promising candidates for ac- 
complishing that, is the artificial version of spin ices re- 
cently produced by Wang et aln. In this system, elon- 
gated magnetic nano-islands are regularly distributed in 
a two-dimensional (2d) square lattice. The longest axis 
of the islands alternate its orientation pointingj-.in the 
direction of the two principal axis of the latticed . The 
magnetocrystalline anisotropy of Permalloy (the mag- 
netic material commonly used to fabricate artificial spin 
ice) is effectively zero, so that the shape anisotropy of 
each island forces its magnetic moment to align along 
the largest axis, making the islands effectively Ising-like. 
Actually, the fabrication and study of this kind of lower 
dimensioH|a4-ianalogues of spin ice have received a lot of 
attention&tJ. Indeed, the ability to manipulate the con- 
stituent degrees of freedom in condensed matter systems 
and their interactions is much important towards advanc- 
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FIG. 1: (Color online) The modified square lattice studied 
in this work. Top: top view of the system. The arrows rep- 
resent the local dipole moments (S a ^) or Sp(i)). Bottom: 
lateral view of the system showing the height offset between 
islands. The original material produced by Wang et aln is 
two-dimensional with h — 0. 



ing the understanding of a variety of natural phenom- 
ena. Particularly in this context, the possibility po£-©b- 
serving magnetic monopoles in artificial spin icesli-fei is 
a timely problem given that these magnetic compounds 
could provide the opportunity to see them up close and 
also watch them move (for example, with the aid of mag- 
netic force microscopy). Very recently, the direct obser- 
vation of these defects in pi artificial kagome lattice was 
reported by Ladak et aZ£3. However, there is a stim- 
ulating challenging for such an observation (or not) in 
artificial square lattices as pointed out in advance. 

In a previous worktil we have pointed out that 
monopoles do not appear as effective low-energy degrees 
of freedom in two-dimensional square spin ices, as they 
do in the three-dimensional materials {Dy,Ho}2Ti20j. 
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FIG. 2: (Color online) Up: In the artificial spin ice proposed 
in Ref. |l0|, the spins obeying the ice rule do not point along 
directions passing by the center of a tetrahedron as they do 
in the natural spin ice compounds. Down: Configurations of 
the spins obeying the ice rule in a tetrahedron in the artifi- 
cial (left) and the natural (right) spin ices. This small dis- 
tortion of the spins configuration causes a residual ordering 
and consequently, an outstanding energetic string connects 
the monopoles in the modified artificial system. 



Due to the antiferromagnetic order in the ground- state, 
the constituents of a pair monopole-antimonopole be- 
come confined by a string which forbids them to move 
independently. However, we have also argued that above 
a critical temperature, the string configurational entropy 
may lose its tension leaving the monopoles free. The 
quantitative analysis of suck-a possible phase transition 
is under current invest igationtS. Meanwhile, other strate- 
gies to find monopoles in synthetic spin ices have been 
proposed. Moller and Moessnertj have suggested a mod- 
ification of the square lattice geometry in which they ar- 
gue that, considering a special condition, the string ten- 
sion vanishes at any temperature. This modification in 
the system produced by Wang et al.u consists of intro- 
ducing a height offset h between islands pointing along 
the two different directionsli3lij (see Fig. |l|; such a sys- 
tem is currently under experimental planningjlj) . Their 
idea comprises basically the following: if h is chosen so 
that the energies of all vertices obeying ice rule become 
degenerate, then, an ice regime is established leaving the 
monopoles "free" to move (indeed, thpe is a Coulom- 
bic interaction between the monopoles)13. For point-like 
dipoles they considered that a degenerate state is ob- 
tained when the interactions between nearest-neighbors 
(Ji) and next-nearest-neighbors {J 2) are equal, leading 
to the following value for the height offset where "free" 
monopoles ope.ur: hi ce ~ 0.419a (where a is the lat- 
tice spacing)L2l. Taking into account the finite exten- 
sion of the dipoles, the height offset diminishes and as 
e = 1 — l/a — >• (I is the length of the island), the 
endpoints of the islands form a tetrahedron, so that at 
h = eaj\/2 the ordering disappears, and the monopoles 



Here we numerically calculate the energetics of the 
ground-states and excitations in the modified square lat- 
tice as a function of h. In our calculations we consider 
point-like dipoles forming the lattice. Although the main 
physical aspects of the system must be correct with this 
approximation, some parameter values (such as magnetic 
charge, string tension, critical height etc) should be quan- 
titatively altered for the realistic case in which I has a 
finite length. On the other hand, since we take into ac- 
count all the long-range dipole-dipole interactions, it is 
expected that our results could better describe the ac- 
tual system. For instance, while in the calculations of 
Refs. |lO| , |T6| , the ground-state changes its configuration 
at h = 0.419a, our results indicate that it occurs at 
h = hi = 0.444a. Besides, we noted that at least one 
of the several configurations that satisfy the ice rule does 
not have the same energy of the "ground-states" (GSi 
and GS2) at this very height, indicating that for h = hi 
the system is not in a completely degenerate state. We 
have also shown that the string tension decreases rapidly 
as h increases but it does not vanish at any value (h < a): 
rather, at h = /ii, its strength reads about 20 times 
smaller than that of the usual case for h = 0. A pos- 
sible cause of the finite strength of the string tension 
even at hi is the fact that, concerning the spin configu- 
rations in a tetrahedron, the artificial spin ice has a slight 
difference with its natural counterpart. For the artificial 
compounds proposed in Ref. the localized magnetic 
moments forming a corner-sharing tetrahedral lattice are 
forced to point along the longest axis of the islands (here, 
x- or ^/-directions, see Fig.||) while in the original 3d spin 
ices, they point along a < 111 > axis (indeed, in this 
case, the magnetic dipoles point along axes that meet 
at the centers of tetrahedra). As a result of this mis- 
match, there is always a single ordered ground state in 
the artificial systems , which is responsible for the resid- 
ual value of the string tension and its anisotropy. An- 
other interesting result obtained here with the point-like 
dipole approximation is that the magnetic charge of the 
monopoles jumps as the system undergoes a transition in 
its ground state. In addition, in general, this strength of 
the interaction between a monopole and its antimonopole 
is anisotropic, depending on the lattice direction and on 
the type of order. However, as expected from the above 
discussions, we note that the system anisotropy dimin- 
ishes as h goes to hi. Actually, as h increases from zero, 
the differences found in the values of the "charges" (as 
distinct directions for the monopoles separation are taken 
into account) decreases, and they tend to disappear as 
h — >• fti, i.e., in the ice regime (nevertheless, h = hi is 
not really an optimal ice regime, at least for point-like 
dipoles). 
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FIG. 3: (Color online) (a) Ground-state configuration for h < 
hi = 0.444a, GSj. Note that this is exactly the same state 
obtained in Refs. [10| , [l4| . (b) Configuration of the ground-state 
GS2 obtained for h > 0.444. In GS2, each vertex has a net 
magnetization but globally the magnetization vanishes. Note 
that the ice rule is manifested in every vertex, (c) Another 
configuration that satisfy the ice rule but has an energy higher 
than the configurations shown in (a) and (b) when h = hi . 



II. THE MODEL AND RESULTS 



We model the system suggested in Refs. yjj,[16| assum- 
ing: the magnetic moment ("spin") of the island is re- 
placed by a point dipole at its center. At each site 
(xi, yi,Zi) of a "square" lattice two spin variables are de- 
fined: S a (j\ with components S x = ±1, S y = 0, S z = 
located at f a = (x{ + a/2, y^ /i), and Sp^ with compo- 
nents S x = 0, S y = ±1, S z = at rp — (xi,yi + a/2, 0). 
Spins pointing along the ^-direction and spins pointing 
along the ^-direction are in different planes, separated 
by a height h (see Fig. El). Hence, in a lattice of volume 
L 2 = n 2 a 2 one gets 2 x n 2 spins (we have studied systems 
with n = 20, 30, 40, 50, 60, 70). Representing the spins of 
the islands by S^, assuming either S a ^ or Spu^ then the 
modified artificial spin ice is described by the following 
Hamiltonian 



Hsi 



Da3 J2 



S% ' Sj 3(S^ • Tij){Sj • Tij) 



, (i) 



where D = fiQfi 2 /4ira 3 is the coupling constant of the 
dipolar interaction. The sum is performed over all 
n 2 (2n 2 — 1) pairs of spins in the lattice for open bound- 
ary conditions (OBC), while for periodic boundary con- 
ditions (PBC) a cut-off radius was introduced when 
> n /2a. 

The results presented here consider a lattice with n = 
70, which contains 9800 dipoles (islands) and PBC. We 
observed exactly the same behavior for OBC and PBC 
and the size dependence of the results is not appreciable. 




FIG. 4: (Color online) The energy per island of the two 
ground-states (GSi and GS2) and of the configuration shown 
in Fig. |3] (c) (in units of D) as a function of h (in units of 
the lattice spacing a). Black circles represent the GSi energy 
while red squares concern GSo and blue diamonds are for the 
configuration shown in Fig. N (c). 



By using a simulating annealing process (see Ref. Q, 
the first thing to notice is that the ground-state configu- 
ration changes for a critical value of h. Indeed, as shown 
in Fig. H (a), for all values h < h\ = 0.444a, the system 
ground-state (hereafter referred to as GS\) has exactly 
the same form as that of the usual case in which h = 0. 
However, for h > h\, the ground-state changes to GS2 
(see Fig. || (b)). Really, as h — >• /ii, the energies of both 
states are comparable, whereas for h > hi the state GS2 
is less energetic (see Fig. |]). Such a result is in qualita- 
tive agreement with findings of Ref. which presents 
the transition at h = hi ce = 0.419a. As expected, both 
configurations obey the ice rule (two spins point in and 
two point out in every vertex) , but while in GS\ the mag- 
netization is zero at each vertex, in GS2 it points diago- 
nally, but with net vanishing magnetization. As shown in 
Fig. [|, the energy of GS\ increases rapidly as h increases 
while the energy of GS2 is constant. Actually, for this lat- 
ter configuration, the horizontal and vertical sub-lattices 
are decoupled. We note that these two ground-states are 
metastable in the sense that they are local minima and 
cannot be continuously deformed one into another with- 
out spending a considerable amount of energy; trying to 
align the dipoles from one state to another costs the in- 
version of two spins by vertex (half of the spins have to 
be inverted in the whole system). This changing has an 
a-dependent energy barrier which is roughly of the or- 
der of 1601} (for h = 0.444a and n = 70), making this 
process much improbable to occur spontaneously. Thus, 
considering the system in the GS\ state and increasing 
continuously the height from h = 0, GS\ may persist 
even for h > hi because of the large energy necessary to 
change to GS2. Besides, in Fig. |]we also present the en- 
ergy of the configuration shown in Fig. || (c), which also 
satisfy the ice rule but has an energy higher than those of 
GSi and GS2 even for h = hi. Consequently, the states 
satisfying the ice rule are not completely degenerate. 

Now, we consider the excitations above the ground- 
state. In the two-in/two-out configuration, the effective 
magnetic charge Q 1 ^ (number of spins pointing inward 
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FIG. 5: (Color online) Three of the four basic shortest strings 
used in the separation process of the magnetic charges. Pic- 
tures (1) and (2) exhibit strings 1 and 2, respectively used 
for h < hi = 0.444a. The red circle is the positive charge 
(north pole) while the blue circle is the negative (south pole). 
For h > hi the ground-state is GS2 and we used a linear 
string-path (not shown above) and a diagonal path (picture 
(3))- 




FIG. 6: (Color online) The monopole "charge" q (see Eq. g) 
obtained analyzing the energy in the separation process of the 
charges for the two string shapes shown in Fig.N for h < hi. 
When h > hi, the charges variation is shown for a linear 
and diagonal string-paths. Here, q is in units of Da while 
h in units of a. Note how the anisotropy of the monopole 
interaction decreases considerably as h —> hi from below. 



minus the number of spins pointing outward on each ver- 
tex is zero everywhere for h < hi (GSi) and for 
h > hi (GS2). The most elementary excited state is 
obtained by inverting a single dipole to generate local- 
ized "dipole magnetic charges". Such an inversion cor- 
responds to two adjacent sites with net magnetic charge 
Q l j^ = ±1, which is_alike a nearest-neighbor monopole- 
antimonopole pairBEj. Following the same method of 
Ref. [l4], it is easy to observe that such "monopoles" can 
be separated from each other without violating the lo- 
cal neutrality by flipping a chain of adjacent spins. We 
choose four different ways they may be separated (see 
Fig. ||). Firstly, using the string shape 1 and starting in 
the ground-state GS\ (for h < h\ = 0.444a) we choose an 
arbitrary site and then the spins marked in dark gray in 
Fig. H are flipped, creating a monopole-antimonopole sep- 
arated by R = 2a. Next, the spins marked in light gray 
are flipped and the separation distance becomes R = 4a, 
and so on. In this case, the string length (X) is related 
to the charges separation distance R by X = AR/2 (the 
monopole and the antimonopole will be found along the 
same horizontal or vertical line). Secondly, we also con- 
sider a string-path of form 2 (for h < hi), making the 
separated monopoles to be found in different lines (di- 
agonally positioned; now we have X = 2R/y/2). More 
two equivalent ways were studied for h > hi in which 
GS2 is the ground-state. In this case, however, dif- 
ferently from the situation in the GS± state, now the 
monopoles can be separated by using a linear string-path 
(so that X = R) without any violation of the ice rule. 
Finally, another monopoles separation studied for GS2 
is the "diagonal-path" (or path 3), in which the charges 
are put in different lines. Our analysis shows that besides 
the Coulombic-type term q(h)/R (where q = f^qiq2 < 
is the a coupling constant which gives the strength of 
the interaction), the total energy cost of a monopole- 
antimonopole pair has an extra contribution behaving 
like b(h)X, brought about by the string-like excitation 
that binds the monopoles, say, 



V(R, h) = q(h)/R + b(h)X(R) + V (h) 



(2) 



where Vo(h) is a a-depended constant related to the 
monopole pair creation (for instance, for h = Vb(0) ~ 
23D and V(a,0) w 29D). The results for the "charge" 
q(h) are shown in Figure ^ for the range < h < a. 
When h < hi, the excitations are considered above GSi, 
and we observe that there is a small h-dependent differ- 
ence in the g- value for paths 1 and 2, which vanishes as 
h — >• hi. At higher heights, h > hi, q is valued with 
respect to GS2, and is a-independent for a linear string- 
path. However, for path 3, it comes back to increase 
as h increases. Therefore, the interaction of a monopole 
with its partner (antimonopole) is anisotropic in artifi- 
cial spin ices. Perhaps, it would be more appropriate to 
redefine things in such a way that q — ^QiQ2a(h,(f)), 
where gig2 = QiQ2&(h,(j)) and the actual value of the 
charges Qi = —Q2 is independent of the angle <p that 
the line connecting the poles makes with the x-axis. In 
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FIG. 7: (Color online) The string tension for the two string 
shapes shown in Fig.^ for ft < hi and for a linear string-path 
and path 3 (diagonal) for ft > hi. The green dot and the 
dashed lines represents an extrapolation of our data. 



this case, the anisotropy of the interaction (coming from 
the background) is implicitly considered in the function 
a (ft, <j)) but its complete expression was not evaluated 
here. Since a(fti,0) tends to be a constant (indepen- 
dent of </>, we set a(fti,0) = 1 and so, only around 
the ice regime (i.e., ft « fti), the interaction tends to 
be isotropic. Thus we can find the genuine strength of 
the magne tic charge in th is artificial compound as being 
Qi = dz \J 47T | q(hi) | J/jLo « ± 1.95/i/a, where we have 
used | g(fti) |= 3.8Da. Just for effect of comparison, us- 
ing some parameters of Ref.^jsuch as a = 320nm, we get a 
charge value which is about 80 times larger than the typ- 
ical value found for the original 3d spin iceslil (or about 
100 times smaller than the Dirac fundamental charge). 
Besides its anisotropy, another interesting fact about the 
Coulombic interaction in the artificial compounds is that 
it jumps at ft = fti. Indeed, at this point, q abruptly 
changes from q < & S.SDa to q > « — 3ADa when the 
linear path is taken into account. Such a discontinu- 
ity may be attributed to the ground-state transition and 
that above GS2 the Coulombic interaction between a pair 
somewhat incorporates the residual magnetization stored 
in each vertex. On the other hand, keeping a diagonal 
separation of the monopoles along the ground state tran- 
sition, the magnetic charge parameter q increases almost 
continuously. 

How the string tension b depends upon ft is shown in 
Fig.^. Note that, while GS\ is the ground-state (ft < fti), 
b diminishes as ft increases. At higher heights, and being 
evaluated over GS2, the tension remains a non- vanishing 
small constant for linear path and turns back to increase 
for diagonal separation (path 3). In general, since b is 
also a function of <\> (i.e.,6(ft, </>)), it is more favorable en- 
ergetically that a pole and its antipole reside at the same 
line in the array. It should be remarked that, near the 
ice regime, 6(fti,0) is almost independent of <j) (almost 
isotropic limit) and its value is reduced around 20 times 
whenever this modified system is compared to its coun- 
terpart at ft = (at zero temperature). In principle, this 
result indicates that free monopoles-do not appear in this 
system. Then, the modified arraytS faces a small obsta- 
cle by the fact that the islands are placed in such a way 



that the spins can not point to the center of a tetrahe- 
dron as they do in the 3d materials. Indeed, as pointed 
out before, the spins in the artificial compound point 
along its edges (see Fig.||); the islands are rigid objects 
that do permit the spins to point only along their longest 
axis. This disparity causes an ordering in the artificial 
material, which diminishes as ft increases; eventually it 
becomes tiny but persists at ft = fti. This persistent or- 
dering contributes for the residual string tension at the 
ice regime and also for the different string tension values 
as the monopoles are located at different angular posi- 
tions in the array. Such a difficulty may be overcome 
when one takes the limit I —> a in the modified array. 



As pointed out in Ref.10, the mechanism responsible for 
the equivalence between the artificial (2d) and natural 
3d spin ices is not operational in d = 2, as it requires 
also the dimensionality of the dipolar interaction to coin- 
cide with that of the underlying lattice. Here, we have a 
d = 3 dipolar (1/r 3 ) and Coulombic ("monopolar", 1/R) 
interactions in a two-dimensional array. Independent of 
this, since the state GS\ is met ast able one could imagine 
if the excitations could be considered to lie in GS\ for ft 
slightly greater than fti. In this case the extrapolation 
of our results indicate that the string tension may vanish 
at ft w 0.502a (see Fig. |). 



III. SUMMARY 

In summary, we have investigated the energetics of the 
modified artificial spin ice expressing several quantities, 
such as ground states energy, magnetic charges and string 
tension, as a function of the height offset ft. Our analysis 
show that the ground-state changes from an ordered anti- 
ferromagnetic to a ferromagnetic one at ft = fti « 0.444a, 
which is in good agreement with the value obtained in 
Refs. |Tc|Jl6| , hi ce « 0.419a. We claim that such a small 
difference comes about from the fact that in these cited 
works, authors assumed equal nearest-neighbor and next- 
nearest-neighbor interactions, whereas we have taken all 
the dipole interactions into account. For the excitations 
above the ground-state we have found that the magnetic 
charges interact through the Coulomb potential added 
by a linear confining term with tension 6(ft), which de- 
creases rapidly as ft increases, from to fti, assuming 
a non- vanishing constant value at higher ft. Actually, 
the system presents an anisotropy that manifests itself in 
both the Coulombic and linear interactions and it tends 
to diminish as ft increases, almost disappearing at ft = fti. 
The source of this anisotropy is a residual ordering, which 
still persists even in the ice regime (at ft = fti for point- 
like dipoles). Ordering and anisotropy may disappear 
completely in the ideal limit I — »• a, ft — >• 0. Another 
interesting result is that the magnetic charge jumps, de- 
pending on the direction in which the monopoles are sep- 
arated, as the system undergoes a transition in its ground 
state. For a separation of the monopoles, vertex by ver- 
tex, along the same line of vertices, which is possible 
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only for the GS2 ground state, the coupling q exhibits 
considerably discontinuity in relation to its limit value 
in the GSi ground state. On the other hand, it tends 
to grow up continuously for the diagonal path along the 
transition. Although the residual ordering leads to a con- 
fining scenario for monopoles, its very small strength, 
whenever h « hi, signalizes a significant tendency of 
monopole-pair unbinding at a critical (optimal) height 
offset, even at zero- temperature. Further improvements 
in model (1), for instance, taking the actual finite-size 
of the dipoles into consideration, could shed some extra 
light to this issue. Additionally, temperature effects may 
also facilitate the conditions for free monopoles. Indeed, 
the string configurational entropy is also proportional^ 
the string size and therefore, at a critical temper at uretll, 



on the order of ba, the monopoles may become free. In 
view of that, for small 6, the monopoles should be found 
unbind at very low temperatures. As a final remark we 
would like to stress that these results show that the back- 
ground configuration of spins has a deep effect in the 
charges interactions, being responsible for the string ten- 
sion, anisotropics and a kind of screening of the charges. 
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